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In this lecture, we give a brief review of what theorists now know, understand, or 
CO ' guess about static and kinetic properties of quark-gluon plasma. A particular attention 

' is payed to the problem of physical observability, i.e. the physical meaningfulness of 

various characteristics of QGP discussed in the literature. 

1 Introduction. 

Ph. 

43 ■ It is well-known for theorists that the physics of QCD medium at high temperatures 

T S> A-qcd differs dramatically from how the system behaves at zero and at low temper- 

^ ■ atures. At low temperatures, we have confinement and the spectrum involves colorless 

C$ ■ hadron states. The lightest states are pions — pseudogoldstone states which appear 

due to spontaneous breaking of (approximate) chiral symmetry of QCD lagrangian. 

At high temperatures, hadrons get "ionized" to quarks and gluons, and, in the 0— 
approximation, the system presents the heat bath of freely propagating colored parti- 
cles. For sure, quarks and gluons interact with each other, but at high temperatures 
the effective coupling constant is small a s (T) <C 1 and the effects due to interaction 
can be taken into account perturbatively 0. This interaction has the long-distance 
Coulomb nature, and the properties of the system are in many respects very similar 
to the properties of the usual non-relativistic plasma involving charged particles with 
weak Coulomb interaction. The only difference is that quarks and gluons carry not the 
electric, but color charge. Hence the name: Quark-Gluon Plasma (QGP). 



lecture at the XXIV ITEP Winter School (Snegiri, February 1996). 

2 We state right now, not to astonish the experts, that there are limits of applicability of perturbation 
theory even at very high temperatures, and we are going to discuss them later on. 
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To describe the properties of QGP is a rather interesting problem of theoretical 
physics lying on cross-roads between the relativistic field theory and condensed matter 
physics. Personally, I had a great fun studying it. Unfortunately, at present only a 
theoretical study of the problem is possible. Hot hadron medium with the temperature 
above phase transition can be produced for tiny fractions of a moment in heavy ion 
collisions but: 

• It is not clear at all whether a real thermal equilibrium is achieved. 

• A hot system created in the collision of heavy nuclei rapidly expands and cools 
down emitting pions. It is not possible to probe the properties of the system 
directly, but only indirectly via the characteristics of the final hadron state. 

• Anyway, the temperature achieved at existing accelerators is not high enough 
for the perturbation theory to work and there are no quantitative theoretical 
predictions with which experimental data can be compared. May be RHIC would 
be better in this respect. 

Thus, at present, there are no experimental tests of nontrivial theoretical predictions 
for QGP properties. The effects observed in experiment such as the famous J/tp - 
suppression (for a recent review see ]|]) just indicate that a hot and dense medium is 
created but says little on whether it is QGP or something else. 

The absence of feedback between theory and experiment is a sad and unfortunate 
reality of our time: generally, what is interesting theoretically is not possible to measure 
and what is possible to measure is not interesting theoretically []. 

Thus, we will not attempt in this lecture to establish relation of the results of the- 
oretical calculations with realistic accelerator experiments. What we will do, however, 
is discussing the relation of theoretical results with gedanken experiments. Suppose, 
we have a thermos bottle with QGP on a laboratory table and are studying it from 
any possible experimental angle. We call a quantity physical if it can in principle be 
measured in such a study and non-physical otherwise. We shall see later that many 
quantities discussed by theorists may be called physical only with serious reservations, 
and some are not physical at all. 



2 Remarks on finite T diagram technique. 

The main point of interest of these notes are the physical phenomena in hot QCD 
system. However, as the main theoretical tool to study them is the perturbation 

3 One of the possible exceptions of the general rule in the field of thermal QCD is a fascinating perspective 
to observe the phenomenon of disoriented chiral condensate at RHIC. But that refers to the temperature 
region T ~ T c ~ 200Mev where the phase transition with restoration of chiral symmetry occurs (for a review 
of physics of phase transition in QCD, see e.g. ||). 
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theory and we want in some cases not only to quote the results, but also to explain 
how they are obtained, we are in a position to describe briefly how the perturbative 
calculations at finite T are performed. 

There are two ways of doing this — in imaginary or in real time. These techniques 
are completely equivalent and which one to use is mainly a matter of taste. Generally, 
however, the Euclidean technique is more handy when one is interested in pure static 
properties of the system (thermodynamic properties and static correlators) where no 
real time dependence is involved. On the other hand, when one is interested in kinetic 
properties (spectrum of collective excitations, transport phenomena, etc.), it is more 
convenient to calculate directly in real time. 

2.1 Euclidean (Matsubara) technique. 

Consider a theory of real scalar field described by the hamiltonian H[<j){x), n(x)]. The 
partition function of this theory at temperature T can be written as 

Z = Tr{e-P H } = [ n#(x)/C[<Mx),<A(x);/3] (2-1) 

where K, is the quantum evolution operator in the imaginary time (5 = 1/T : 

/C[<//(x),</>(x);/?] =5>«[^(*)]*n[#x)]e-^ (2.2) 

n 

\I/ n are the eigenstates of the hamiltonian. One can express the integral in RHS of 
Eq.(|2.1|) as an Euclidean path integral: 

Z = J JJd0(x,r)exp j-jj^dr J dx£[0(x,r]| (2.3) 

where the periodic boundary conditions are imposed 

<Xx,t + /?) = 0(x,t) (2-4) 
A thermal average < O >t of any operator 0[(f)(x, r)] has the form Q 

< O > T = Z- 1 J ]Jd(f)(x,T)0[(i)(x,T)}exp ^ - J%t J dx£[0(x,r]| (2.5) 

One can develop now the diagram technique in a usual way. The only difference 
with the zero temperature case is that the Euclidean frequencies of the field </>(x, r) are 
now quantized due to periodic boundary conditions 



p% = ImnT (2.6) 
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with integer n. To calculate something, one should draw the same graphs as at zero 
temperature and go over into Euclidean space where the integrals over Euclidean fre- 
quencies are substituted by sums: 

/(gr/W-^E/l^/^p) (2.7) 

The same recipe holds in any theory involving bosonic fields. In theories with fermions, 
one should impose antiperiodic boundary conditions on the fermion fields ip(x,r) (see 
e.g. H for detailed pedagogical explanations), and the frequencies are quantized to 

p% = m{2n + 1)T (2.8) 

An important heuristic remark is that, when the temperature is very high, in many 
cases only the bosons with zero Matsubara frequencies po = contribute in < O >t- 
The contribution of higher Matsubara frequencies and also the contribution of fermions 
in < O >t become irrelevant and, effectively, we are dealing with a 3- dimensional 
theory. As was just mentioned, it is true in many, but not in all cases. For example, 
it makes no sense to neglect fermion fields when one is interested in the properties of 
collective excitations with fermion quantum numbers. For any particular problem of 
interest a special study is required. 



2.2 Real time (Keldysh) technique. 

Matsubara technique is well suited to find thermal averages of static operators 0(x). 
If we are interested in a time-dependent quantity, there are two options: i) To find 
first the thermal average < 0(x, r) > for Euclidean r and perform then an analytic 
continuation onto the real time axis. It is possible, but quite often rather cumbersome, 
ii) To work in the real time right from the beginning. The corresponding technique 
was first developed in little known papers || and independently by Keldysh [J?J who 
applied it to condensed matter problems with a particular emphasis on the systems out 
of thermal equilibrium. It was fully apprehended by experts in relativistic field theory 
only in the beginning of nineties. 

As we do not really need the full-scale Keldysh technique in what follows, we 
will not discuss it here, referring a reader to the textbook ||, to a hard-to-read, but 
extensive review ||, and to our papers |]i~0| , |TT| , |i~2|| where the real time technique was 
applied for studying the properties of QGP. We only mention here that at finite T one 
should distinguish carefully between the retarded, advanced and mixed components of 
the Green's functions which can be "organized" in a 2 x 2 matrix. 

It suffices often to use a simplified version of the real time technique due to Dolan 
and Jackiw [13]. Let us find the tree propagator of real scalar field at finite temperature: 

< 0(x)0(O) > T = e ~ PEn < n\<p(x)(j)(0)\n > (2.9) 
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Introduce as usual a finite spatial volume V and decompose 

<t>( X ) = Y -pL= 

where ap~ and a p are the creation and annihilation operators. At zero temperature 
only the vacuum average contributes and we can use the fact a p |0 >= to obtain a 
usual expression for the propagator. At finite temperature, the excited states \n > 
with a p | re contribute in the sum ( |2.9[ ), and an additional contribution to the 
propagator arises. We obtain 

/ < 0(x)0(O) > T e*** = + 2vr5(p 2 - m 2 )n B (\p \) (2.11) 

J p z — m A + lO 

where 

n B (e) = (2.12) 

is the Bose distribution function. An additional thermal contribution reflects the pres- 
ence of real bosons on mass shell in the heat bath. A fermion propagator is derived in 
a similar way and an additional contribution is proportional to the Fermi distribution 
function np(e) = \ j(eP e + 1). 

The recipe is to substitute the Dolan-Jackiw propagators (and in cases when it does 
not work — the full-scale matrix Keldysh propagators) in the loop integrals for the 
usual Feynman graphs. The integrals are done now over real frequencies. 



a p e-^ +ipx + a p e 



iuit— ipx 



(2.10) 



3 Static Properties of QGP: a Bird Eye's View. 
3.1 Thermodynamics 

The basic thermodynamic characteristic of a finite T system is its free energy. We have 
in the lowest order for a pure gluon system 



F 9 



-ThxZ 



-Tin 



EE 



-Mp\ 



2(JV c a -l)" 



(3-1) 



\n=0 



where 2(iV ( ? — 1) is the number of degrees of freedom of the gluon field (the factor 2 
comes due to two polarizations) and the sum J2 n e _,fln ' p ' is the free energy of a single 
boson field oscillator with the frequency u = |p|. Trading the sum for the integral, we 
obtain for the volume density of the free energy 

d 3 p 



F9 
~V 



2(N* - 1)T 



(27T) 



In 



-0|P| 



7T 2 T 
'^5 



W ~ 1) 



(3-2) 
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This is nothing else as the Stefan-Boltzman formula multiplied by the color factor 
N% — 1. The quark contribution is obtained quite similarly: 



F q = -ThxZ = -Tin 



T-r / I \ 4N c Nf 

TJ(l + e -^|p| x 



p 



(3.3) 



where 4N c Nf is the number of degrees of freedom and the Pauli principle is taken into 
account. We obtain 



V f J (2tt) 3 L J 180 



N c N f (3.4) 



All other thermodynamic quantities of interest can be derived from the free energy by 
standard thermodynamic relations. For example, the pressure just coincides with the 
free energy with the sign reversed. The energy density is 

dF 

E = F - T 9T (3 ' 5) 
For massless particles in the lowest order the relation 

E = 3P = -3F (3.6) 

holds. 

One of the simple and instructive exercises which can be now done is looking at 
the limit N c — » oo. We see that the energy density becomes infinite in this limit. That 
means that if we start to heat the system from T = 0, we just cannot reach the QGP 
state — to this end, an infinite energy should be supplied ! 

This physical conclusion can also be reached if looking at the problem from the low 
temperature end. The spectrum of QCD in the limit N c — > oo involves infinitely many 
narrow states. The density of states grows exponentially with energy [] 

p{E) oc e cE (3.7) 

That means that the partition function 

Z ~ J P {E)e~ E/T dE (3.8) 



4 One of the way to see it is to use the string model for the hadron spectrum. A string state with large 
mass is highly degenerate. The number of states with a given mass depends on the number of ways p(N) the 
large integer N ~ M 2 /a (a is the string tension) can be decomposed in the sums of the form N — J^i n i ( see 
e.g. [0). p(N) grows exponentially with N. That does not mean, of course, that in real QCD with large 
N c the spectrum would be also degenerate. Numerical calculations in QCD2 with adjoint matter fields show 
that there is no trace of degeneracy and the spectrum displays a stochastic behavior [|l5) . And that means 
in particular that there is little hope to describe quantitatively the QCD spectrum in in the limit N c — > 00 
in the string model framework. But a qualitative feature that the density of states grows exponentially as 
the energy increases is common for the large N c QCD and for the string model. 
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is just not denned at T > c . There is a Hagedorn temperature Tjj = c above 
which a system cannot be heated [EfJ. 

When N c is large but finite, no limiting temperature exists (It is seen also from 
the low temperature viewpoint: at finite N c the states have finite width and starting 
from some energy begin to overlap and cannot be treated as independent degrees of 
freedom), and one can bring the system to the QGP state when supplying enough 
energy. But when iV c is large, the required energy is also large. That suggests (though 
does not prove, of course) that at large finite N c a first order phase transition with 
a considerable latent heat takes place. This conjecture is supported by the lattice 
measurements, but we will not pursue this discussion further and refer an interested 
reader to our review ||. 



3.2 Debye screening. 

Consider the gluon polarization operator in QGP with account of thermal loop cor- 
rections. It is transverse, kfjl^ u (k) = 0. At zero temperature, transversality and 
Lorentz-invariance dictates the form Hn V (k) = Il(k 2 )(g fiu — k^k v /k 2 ). At finite T, 
Lorentz-invariance is lost and the polarization operator presents a combination of two 
different (transverse and longitudinal) tensor structures. Generally, one can write 



n o = ii^,|k| 

k{U) . 
kT 



n « = TTj2 n i( w ) l k 



U i:j = -U t (u, |k|)(<% - k t k,/\k\ 2 ) + ^^n,( w , |k|) (3.9) 

Consider first the longitudinal part of the polarization operator in the kinematic region 
where u is set to zero in the first place after which k = |k| is also sent to zero. By the 
reasons which will be shortly seen, we denote this quantity m 2 D : 

m% = ]imUi(0,k) (3.10) 

To understand the physical meaning of this quantity, consider the correlator 

< ^g(x)4(0) >~ 5 ab f - 9 dk 9 e* kx oc e - mD ^ (3.11) 
J k z + m D 

Thus, mo coincides with the inverse screening length of chromoelectric potential Aq. 

There is a clear analog with the usual plasma. A static electric charge immersed 
in the plasma is screened by the cloud of ions and electrons so that the potential falls 
down exponentially oc exp{— r/ro} where the Debye radius is given (in an ordinary 
non-relativistic plasma) by the expression (see e.g. ||]) 
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a) b) c) d) 

Figure 1: Gluon polarization operator in one loop. 



where n is the electron and ion density, Z is the ion charge (so that the second term 
describes the ion contribution), and, to avoid unnecessary complications, we assumed 
that the electron and ion components of the plasma have the same temperature T 
Calculating the one-loop thermal contribution to the gluon polarization operator (see 
Fig.l), one can easily obtain an analogous formula for QGP: 

ml = r» 2 = 9 ^f (N c + ^) (3.13) 

where the first term describes the screening due to thermal gluons and the second 
term — the screening due to thermal quarks. The result ( |3.13 ) was first obtained by 



Shuryak |17|| . It has the same structure as Q3,12|) (Note that the density of particles 
in ultrarelativistic plasma is expressed via the temperature, n oc T 3 ). It is worth 
mentioning that the Debye screening is essentially a classical effect and not only quarks 
but also gluons result in screening rather than antiscreening. That should be confronted 
with the famous antiscreening of the charge in Yang-Mills theory at zero temperature 
due to quantum effects. 

In QED, the correlator < Aq(x.)Aq(0) >t is a gauge invariant object and also 
the notion of charge screening is unambiguous and well defined — there are classical 
electric charges and one can measure the potential of such a charge immersed in plasma 
by standard classical devices. Not so in QCD. We do not have at our disposal classical 
color charges due to confinement and the experiment of measuring the chromoelectric 
field of a test color charge cannot be carried out even in principle. Also the gluon 
polarization operator which enters the definition ( |3,11| ) of a Debye screening mass is 
generally speaking gauge-dependent. Thus, a question of whether the Debye screening 
mass is a physical notion in a non-abelian theory is fully legitimate. 

The answer to this question is positive. But with reservations. 

Note first of all that though the gluon polarization operator is generally gauge- 
dependent, the Debye screening mass defined in ( p.lO| ) does not depend on the gauge 
in the leading order (we shall discuss also non- leading corrections in the next section). 



S 



That suggests that the result ( 3.13j ) has an invariant physical meaning. 
Indeed, one can consider a gauge-invariant correlator 



C(x) =< P(x)P*(0) > T (3.14) 

where 

P(x) = irTre X p|i^Ao(x,r)dr| (3.15) 

is a gauge-invariant operator (the so called Polyakov line, a special thermal version of 

the Wilson line |l8|]). If r = |x| is not too large (the exact meaning of this will be 
specified later), the correlator has an exponential behavior 

C(r) oc [G 00 (r)] 2 oc exp{-2m D r} (3.16) 



The correlator ( 3.14j ) can be attributed a physical meaning. — TTnC(r) coincides with 



the change of free energy of the system when putting there a pair of heavy quark and 
antiquark at distance r, the averaging over color spin orientations being assumed Jl9|]. 

Also we shall see in the next section that the perturbative corrections ~ g 3 T 4 to a 
perfectly well defined and physical quantity — the free energy — depend directly on 
the value of mo- 

Thus, the notion of Debye screening mass in QGP is physical, though, unfortu- 
nately, not to the same extent as it is in the usual plasma. The correlator ( 3.14j ) with 



the correlation length of fractions of a fermi cannot be directly measured even in a 
gedanken laboratory experiment — at least, we cannot contemplate such an experi- 
ment. But it can be measured in lattice numerical experiments which is almost as 
good. We shall see later that a number of other characteristics of QGP have a similar 
semi-physical status. 



3.3 Magnetic screening. 

Debye mass describes the screening of static electric fields. In abelian plasma, magnetic 
fields are not screened whatsoever. Such a screening could be provided by magnetic 
monopoles, but they are not abundant in Nature. The absence of monopoles is tech- 
nically related to one of the Maxwell equations d{Bi = 0. In the non-abelian case, the 
corresponding equation reads T>{Bi = where T> is a covariant derivative. Thus, gluon 
field configurations with local color magnetic charge density p m ~ d{Bi (with usual 
derivative) are admissible. The presence of such configurations in the gluon heat bath 
results in screening of chromomagnetic fields. 

Let us see how it comes out in a perturbative calculation. Consider the one-loop 
graph in Fig. la which contributes to the polarization operator of the spatial compo- 
nents of the gluon field Ily (0, k) oc IT(0, k). At high temperatures, it suffices to take 
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into account only the lowest Matsubara frequency (with u = 0) in gluon propagators. 
We have 

n t (0, k) ~ g 2 k 2 T J p2( ^ k)2 ~ g 2 Tk (3.17) 

The numerical coefficient can be explicitly calculated, but it depends on the gauge and 
makes as such a little sense. The loop integral is determined by the low momentum 
region (pjchar ~ k <C T (as we are interested in the large distance behavior of the 
gluon Green's function, we keep k small). Note that the Feynman integral for 11/(0, k) 
has a completely different behavior being saturated (in the leading order) by the loop 
momenta |p| ~ T — it is a so called hard thermal loop. 
The transverse part of the gluon Green's function is 

a ' M ~ vrhoj) (3 - 18) 

We see that at k ~ g 2 T the one-loop contribution to the polarization operator is of the 
same order as the tree term k 2 . One can estimate a two- loop contribution to 11^(0, k) 
which is of order (g 2 T) 2 . The factor T 2 here comes from two loops (we use the rule 
( ]2.7|) and take into account only the contribution of the lowest Matsubara frequency 
in each loop). Similarly, a three-loop contribution to 11^(0, A;) is of order (g 2 T) s /k, 
a four-loop contribution is of order (g 2 T) 4 /k 2 etc. (The growing powers of k in the 
denominator are provided by infrared 3-dimensional loop integrals. Infrared integrals 
may also provide for a logarithmic singularity in external momentum k in the two- 
loop contribution which is of no concern for us here). At k ~ g 2 T all contributions are 
of the same order and the perturbation theory breaks down [ |20f . 

There is an alternative way to see it. Let us write down the expression for the 
partition function of the pure glue system at finite temperature: 

Z = J n^(x,T)exp{-^ J*dT J d 3 x(F^ (3.19) 

When T is large and (3 small, the Euclidean time dependence of the fields may be 
disregarded. Also the effects due to Ao(x,r) can be disregarded — time components 
of the gluon field acquire the large mass mc ~ gT S> g 2 T and decouple. We are left 
with the expression 

Z = J n^?(x)exp |-_L J d 3 x(i^.) 2 j (3.20) 

A theory with quarks is also reduced to ( 3.20 ) in this limit — the fermions have high 
Matsubara frequencies ~ T and decouple. The partition function ( |3.20| ) describes a 
non-linear 3D theory with the dimensional coupling constant g 2 ~ g 2 T PL No pertur- 



5 One should be careful here. It would be wrong to use the expression (3.2C) for calculating, say, the 
free energy of QGP at large T. The latter takes the contributions from hard thermal loops (involving 
also fermions!) with momenta of order T which are not taken into account in Eq.(3.2C). ( 3.20| ) should be 
understood as an effective theory describing soft modes with momenta of order g 2 T . 
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bative calculations in this theory are possible. 

In QED, there is no magnetic screening (and the effective 3-dimensional theory is 
trivial). Presumably, perturbative calculations in hot QED can be carried out at any 
order, though this question is not absolutely clear. 

What is the physical meaning of this non-abelian screening ? Can it actually be 
measured ? 

The way we derived it, the magnetic screening shows up in the large distance 
behavior of the spatial gluon propagator Gt(r) oc exp{ik*r} where k* is the solution 
to the equation 

k 2 +n t (0,fc) = (3.21) 

There is no reason to expect that k = is a solution — the behavior of IT(0, k) in the 
limit k — > where the pertubation theory does not work is not known, but one can 
tentatively guess that it tends to a constant ~ (g 2 T) 2 (As we have seen, g 2 T is the 
only relevant scale in this limit). If so, k* ~ ig 2 T and Gt(r) falls down exponentially 
oc exp{— Cg 2 Tr} at large distances. Hence the term "magnetic screening". 

Unfortunately, the gluon polarization operator is a gauge-dependent quantity and, 
would the God provide us with the exact expression for Yl t (0,k) in any gauge, even 
He could not guarantee that the solution of the dispersive equation ( |3.21| ) would be 
gauge-independent []. 

What one can do, however, is to consider the correlator of chromomagnetic fields 
Caft(x) = < B a (x)B 6 (0) >t- This correlator also depends on the gauge in a non- 
abelian theory, but the gauge dependence amounts to rotation in color space: C a b(x) — > 
^a'a( x )C a ' b '(x)f2f/f>(0) an d cannot affect the exponential behavior of the propagator. 
This is the way the magnetic mass is usually measured on the lattices. 

One can do even better considering gauge- invariant correlators, the simplest one is 
< G 2 u (x)G 2 u (0) > T - Then the quantity 

M = - P l™ I ln < G U*) G U°) >T ( 3 - 22 ) 

provides an invariant definition of the magnetic screening mass. The indefinite article is 
crucial here. Choosing other gauge-invariant correlators, one would get other invariant 
definitions. For example (and this will be important in the following discussion) , 



the true correlation length of the correlator of Polyakov loops ( 3.14 ) is also of order 



(g 2 T)~ l . The asymptotics (|1|) is an intermediate one and holds only in the region 

(gT)- 1 « r « (g 2 T)~ 1 (3.23) 

What was wrong with our previous derivation ? The matter is we took into account 
earlier only the thermal corrections to the gluon propagator and tacitly neglected cor- 
rections to the vertices. Ab accurate analysis plj shows that this is justified in the 



3 We shall return to the discussion of this important point in the following sections. 
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0=0 



Figure 2: 



range ( |3.23| ) but not beyond. An example of the graph providing the leading contribu- 
tion in < P(x)P*(0) > T at r » (g 2 ^ 1 is given in Fig.2. [] 

Thus, the "experimental status" of the magnetic mass Q3.22 ) and of other similar 



quantities is roughly the same and even better than for the Debye mass. The expo- 
nential behavior cx exp{— fir} of the correlators is expected to hold at any large r and 
the value \i can in principle be determined in lattice experiments with any desired 
accuracy (now the accuracy is very poor due to a finite size of available lattices, but it 
is not a question under discussion here). On the other hand, the Debye screening mass 
cannot be determined with an arbitrary accuracy due to the finite range ( fTgj) where 



the asymptotics ( 13,16 ) holds (being sophisticated enough, an invariant definition of 



Debye mass still can be suggested for N c > 3 — see the discussion in the following sec- 
tion). And, as was also the case for the Debye mass, we cannot invent any laboratory 
experiment where the magnetic mass could be directly measured. 

4 Static Properties of QGP: Perturbative Cor- 
rections. 



4.1 Debye mass. 

To provide a smooth continuation of the discussion started at the end of the previous 
section, we consider first higher-loop effects in the Debye mass. Note first of all that 

7 The asymtotics of the correlator of Polyakov loops at large enough r is determined by the magnetic 
photon exchange also in abelian plasma. Two magnetic photons can be coupled to P(x) via a fermion loop. 
In abelian case, magnetic photons are massless and, as a result, the correlator has a power asymptotics 
cx 1/V 6 [§2| . Physically, it corresponds to Van-der-Vaals repulsion between the clouds of virtual electrons 
and positrons formed near two heavy probe charges. Seemingly, the asymptotics oc 1/r 6 for the Polyakov 
lines correlator in QGP found in ||^] has this origin. But in non-abelian case, the asymptotics becomes 
exponential when taking into account the magnetic screening effects. 
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the definition ( 3.1C| ) is not suitable anymore. Higher-loop corrections Amc ~ g 2 T to 
the Debye mass as defined in Eq.( 3.10| ) depend on the gauge \24\- A better way is to 
define the Debye mass as the solution of the dispersive equation k 2 + II; (0, k) = 0. In 
other words, we define [p5|] 



m 2 D = lim 11/(0, k) (4.1) 



fc 2 ^— m' 2 D 



The longitudinal part of the gluon Green's function 

1 



Gi(0,k) 



fc 2 + IL,(0,A;) 

has then the pole at k = imrj. It is conceivable that the pole position is gauge- invariant 
even though the Green's function itself is not. (See, however, a discussion in the 
following section. There are cases when formal arguments displaying gauge-invariance 
of the pole position fail in higher orders of perturbation theory due to severe infrared 
singilarities. It is not clear for us whether they really work in the case of Debye mass.) 
In the next-to-leading order, the (gauge-invariant) result is [25] 



Am 2 n N P I 3 , 1 



'■D 



g In- (4.2) 



m 2 D 2tt y N c + Nf/2 g 

We see that the correction is non-analytic in coupling constant. The non-analyticity 

appears due to bad infrared behavior of the loop integrals — they involve a logarithmic 

infrared divergence and depend on the low momenta cutoff which is of order of magnetic 

mass scale g 2 T. Thus, we have 

2 [™d dp m D 1 
Am D oc / — = In ~ in — 

J Umag P ^mag 5 

The magnetic infrared cutoff is actually provided by higher-order graphs — the or- 
ders of perturbation theory are mixed up and a pure two-loop calculation is not self- 
consistent. As we have seen, [i ma g cannon be determined analytically which means 
that the correction ~ g 2 T without the logarithmic factor in the Debye mass cannot be 
determined analytically. 

Also the correction ~ g 2 T cannot be "experimentally observed". Really, we have 
seen that the invariant physical definition of refers to the correlator of Polyakov 
loops ( 3.14j ). The latter displays the exponential behavior ( |3.16 ) in the limited range 
( |3.23 ). But that is tantamount to saying that the correction ~ g 2 T in the Debye mass 
cannot be determined from the correlator ( |3.14[ ). To do this, one should probe the 
distances r 3> {g 2 T)~ l where the correlator has a completely different behavior being 
determined by the magnetic scale. The correction (|4.2| ) is still observable, however, 
due to a logarithmic enhancement factor. 

To be more precise, the correlator in the range (|3.23|) has the form pl| 



< P(x)P*(0) > T oc exp \ -2m$r 




2tt V N c + Nf/2 



g ln(m^V) 



(4.3) 
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where is the lowest order Debye mass ( |3~T3|) . At finite r <C (g 2 T) , the infrared 



cutoff in the loop integrals is provided by r 1 rather than \i ma g and we have 

In mi) > ln(m£)r) 

^mag 

The Fourier image of the correlator ([D|) does not have a singularity at finite k what- 
soever. Thus, the Debye mass pole in the gluon propagator does not show up as a pole 
in the gauge- invariant correlator of Polyakov loops. 

But, anyway, in the theoretical limit g(T) — > there is a range of r where the 



condition r <C {g 2 T) 1 is fulfilled so that the intermediate asymptotic law (O) holds 



whereas the correction in the exponent oc grn^r ln(m£>r) is large compared to 1 and 
can be singled out in numerical lattice experiment. 

There are two interesting recent proposals to define Debye mass non-perturbatively 
in a gauge- invariant way |p2| . The first one is to consider the correlator of imaginary 
parts of Polyakov lines in a complex representation which is invariant under the action 
of the center of the group Zjy (say, the decouplet representation in SU(3)) The 
point is that P* — Pio is odd under Euclidean time reversion and cannot be coupled 
to magnetic gluons. That means that the correlator 

<P 1 * (x)-P 10 (x), P{ (P)-P 10 (P)>t (4.4) 

exhibits the Debye screening falloff oc exp{— 3m£>r} even at arbitrary large r. This 
definition does not work, however, for SU (2) where the Polyakov line is real in any 
representation. 

Another suggestion was to study the behavior of large spatial Wilson loops in 
adjoint color representation. At large distances, the theory is effectively reduced to the 
3- dimensional YM theory ( p. 20 ). Adjoint color charges in this theory are not confined 



but rather screened, and the Wilson loop exhibits the perimeter law behavior 

W{C) oc exp{— m* x perimeter(C)} (4.5) 

One can identify m* (which is of order g 2 T) with a non-perturbative correction to the 
Debye mass. The problem here that m* has no trace of the lowest order contribution 
( |3.13 ). It involves a part of the perturbative correction ( |4.2[) ~ g 2 Tln | due to gluon 



loops but not a similar contribution from fermion loops. Still m* certainly has an 
invariant meaning and is as such an interesting quantity to study. 



8 The authors of Ref. |22| argued the necessity to consider only the Zm - invariant representations saying 
that otherwise the correlators give zero after averaging over different "Zn - phases" . Actually, such phases 
do not exist in Nature and a proper agrumentation should be that the correlators which are not invariant 
under Z/v transformations just do not have a physical meaning pq|. 
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Figure 3: Free energy in two loops. 




Figure 4: An infrared- divergent 3-loop graph in free energy. 



4.2 Free energy. 

This is the most basic and physical quantity of all. May be this is the reason why 
perturbative corrections are known here with record precision. 

The correction ~ g 2 T has been found by Shuryak [17]. To this end, one should 



calculate the two-loop graphs depicted in Fig. 3. The behavior of Feynman integrals in 
this order is quite benign and no particular problem arises. 

However, on the threedoop level, the problem crops up. The contribution of the 
graph depicted in Fig. 4 is infrared divergent: 



AF Fig 



(4.6) 



That means that a pure 3doop calculation is not self- consistent and, to get a finite 
answer, one has to resum a set of infrared-divergent graphs in all orders of pertur- 
bation theory. This is, however, not a hopeless problem, and it has been solved by 
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Figure 5: Ring diagrams. Grey circles stand for a set of one-loop graphs in Fig.l. 



Kapusta back in 1979 [27]. The leading infrared singularity is due to the so called 
"ring diagrams" depicted in Fig. 5. 

The sum of the whole set of ring diagrams has the form 



pnng ^ T / d 3 fc In 



1 + 



n f (o,fc) 
k 2 



ih(Q,fc) 
k 2 



(4.7) 



The expansion of the integrand in 11/(0, k) /k 2 restores the original infrared-divergent 
integrals. However, the whole integral is convergent being saturated by momenta of 
order k ~ [11; (0, 0)] 1 / 2 ~ gT. Neglecting the fc-dependence in the polarization operator, 
we obtain 



Fa 



d 3 k 



D 

k 2 



Tm 



D 



3^4 



(4.8) 



We see that the correction is non-analytic in the coupling constant a s which exactly 
reflects the fact that the individual graphs diverge and the orders of perturbation theory 
are mixed up. Note, however, that infrared divergences here are of comparatively 
benign variety — the integral depends on the scale k char ~ gT and we are far from the 
land of no return k ~ g 2 T. That is why an analytic determination of the coefficient in 
Q4.S{ ) is possible. 

The next correction has the order ~ g 4 T 4 ln(l/g). It comes from the same ring 
graphs of Fig. 5 where now the term ~ g 2 Tk in 11/(0, k) should be taken into account. 
This correction was determined by Toimela [p8| . 

Presently, the correction ~ g 4 T 4 (without the logarithmic factor) |2J| and the 
correction ~ g 5 T 4 pi p| are known. This is the absolute limit beyond which no 
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perturbative calculation is possible — similar ring graphs as in Fig. 5 but with magnetic 
gluons would give the contribution oc Tp? mag ~ g 6 T 4 in the free energy. And, as we 
have seen, fi mag cannot be determined perturbatively. 

Collecting all the terms, the following result is obtained 31] 

a s (/^ 3/2 



where 



8vr 2 T 4 



45 



F + F 2 -^- + F 3 



TT 



TT 



+ 



TT 



a 



F 4 — +F 5 — +0(ailna 



5/2 



TT 



21 



32 



(4.9) 



(4.10) 



15 
T 



1 + — JV> 
12 ' 



(4.11) 



iVf \ 3/2 

••lo i : i + * 



(4.12) 



F 4 = 237.2 + 15.97JV/ 



135 



0.413AT) + ^ h + 12 



165 



1 + — JV> 
12 ; 



In 



vr V 6 
33 ' 



In 



2vrT 



(4.13) 



1 + 



6 



1/2 



-799.2 - 21.96iV 



+- 



485 



N f 

1 + f 



1 - —N f ) In 



1.926JVJ 



(4.14) 



2 V 6 J \ 33 V 2vrT 

The expressions (|4.10| ) - ( |4. 14j) are written for N c = 3. The coefficients like 237.2 are 
not a result of numerical integration but are expressed via certain special functions. 

A nice feature of the result ( [4.9D is its renorm-invariance. The coefficients F4 and 
F$ involve a logarithmic /U-dependence in such a way that the whole sum does not 
depend on the renormalization scale fi. 

Let us choose [i = IttT (this is a natural choice, 2irT being the lowest non-zero 
gluon Matsubara frequency ). In that case, we have 



1 - 0.9q s + 3.3a 3 / 2 + (7.1 + 3.5 In a s )a 2 s - 20.8a 5 J 2 



(4.15) 



5/2 - 

Note a large numerical coefficient at a s ■ It is rather troublesome because the cor- 

5/2 

rection ~ a s overshoots all previous terms up to very high temperatures and, at 
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temperatures which can be realistically ever reached at accelerators, makes the whole 
perturbative approach problematic. 



Take T ~ 0.5 GeV (this is the temperature one can hope to achieve at RHIC 32]). 
Then 2irT ~ 3 GeV and a s ~ 0.2. (We use a conservative estimate for a s following 
from T physics |33|]. Recent measurements at LEP favor even larger values.). The 
series ( |4 .15 ) takes the form 



F = F [l - .18 + .3 + .06 - .37 + . . .] (4.16) 

5/2 

which is rather unsatisfactory. We emphasize that the coefficient of a s is rather 
trustworthy being obtained independently by two different groups. 

The last remark is technical. The result ( |4.9| ) was obtained in Euclidean technique. 
There is a real time calculation which correctly reproduces the two-loop term ~ g 2 T 4 
in free energy ||, but nobody so far succeeded in calculating in this way the terms 
~ g 3 T 4 and higher. Certainly, real-time technique is not very suitable for calculation 
of static quantities, and one way to get the result is good enough, but, to my mind, it 
is an interesting methodical problem. 



5 Collective excitations. 

One of striking and distinct physical phenomena characteristic of usual plasma is a 
non-trivial dispersive behavior of electromagnetic waves. In contrast to the vacuum 
case where only transverse photons with the dispersive law u = |k| propagate, two 
different branches with different non-trivial dispersive laws to±(k) and uj^(k) appear in 
plasma. The value o;_i_(0) = ^|i(0) = oj p i characterizes the eigenfrequency of spatially 
homogeneous charge density oscillations and is called the plasma frequency. 

A similar phenomenon exists also in QGP. The spectrum of QGP involves collective 
excitations with quantum number of quarks and gluons. Like in usual plasma, there are 
transverse and longitudinal branches of gluon collective excitations (alias, transverse 
and longitudinal plasmons) and their properties on the one-loop level are very similar 
to the properties of photon collective excitations in usual plasma. A novel feature is the 
appearance of non-trivial fermion collective excitations (plasminos). But, again, they 
are not specific for a non-abelian theory and appear also in ultrarelativistic e + e~ plasma 
(in the limit when the electron mass can be neglected compared to the temperature- 
induced gap in the electron spectrum cx eT). 

5.1 One-loop calculations. 

The dispersive laws of quark and gluon collective excitations can be obtained via solu- 
tion of a non-abelian analog of the Vlasov system involving the classical field equations 
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in the medium and Boltzmann kinetic equation 34]. This way of derivation makes 
analogies with usual plasma (where the Vlasov system is a standard technique) the 
most transparent. 

I will outline here another way of derivation which is more conventional and more 
easy to understand for a field theorist. This is actually the way the results were 
originally derived |}(], 37]. 

Consider the gluon Green's function in a thermal medium. As was mentioned in 
sect. 2, at T ^ 0, different kinds of Green's function exist. To be precise, we are 
considering now the retarded Green's function 



G>(s) 



R 



l <e(t)[A a Jx),A b l/ (0)] > T 



(5.1) 



which describes a response of the system on a small perturbation applied at t = at 
some later time t > 0. The Fourier image of ( |5.1| ) is free of singularities in the 
upper uj half-plane. The poles of GQ,(w,\s) correspond to eigenmodes of the system 
and exactly give us the desired spectrum of gluon collective excitations. The dispersive 
equation 



splits up in two: 



det||G* (w 







k 2 -IL t (u,k) = 
k 2 + Ui(u>, k) = 



(5.2) 



(5.3) 



with Ui t t(u), k) being defined in Eq. (|3.9D . The solutions to the equations (|5.3|) give two 
branches of the spectrum. 

The explicit one-loop expressions for Ui(uj,k) and Ut(uj,k) obtained by the calcu- 
lation of the graphs in Fig.l in the limit u,k <C T are |3^, [37j 



Tl l {uj,k)=Zuj 2 [l-F{uj/k)\ 



U t (u,k) = -u pl 



+ 



' J P i 
k 2 -u 2 



k 2 



F(u/k) 



(5.4) 



where 



3 V 2 



(5.5) 



9 To find quark dispersive laws, one has to write down generalized kinetic equations for spinor densities. 
Such equations were first introduced in |}5| when studying the problem of goldstino dispersion law in a 
supersymmetric thermal medium. 

10 As in commonly used gauges G ab cx S ab , we shall suppress color indices in the following. We have 
retained them here just to make clear that we are dealing with a commutator of Heisenberg field operators, 
not with a commutator of classical color fields. 
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Figure 6: Plasmon spectrum in one loop, 
is the plasma frequency and 

F(x) = |ln|±j (5.6) 

The behavior of dispersive curves is schematically shown in Fig.6. At k = 0, 
^±(0) = wm(O) = u; p ;. Then the two branches diverge: 



« gT) = u 2 pl + -k 2 
3 



u, 2 (fc« 5 T)=^ + ^ 2 (5.7) 



At k ^> gT both branches tend to the vacuum dispersive law 

uj\{k > 5 T) = k 2 + 3 

w 2 (fc > #T) = /c 2 



2 

2/c_ 2 

3a; 2 



1.+ l(<xp(-T^- 2 l ( r '- 8 ) 

We see that 0J\\{k) approaches the line u = k exponentially fast. 

The dispersive laws for quark collective excitations are obtained from a similar 
analysis of the quark Green's function. At T ^ 0, the fermion polarization operator 
involves two tensor structures 0^7° — k7 and W7 + k7 which gives rise to two dispersive 
branches. We will call the branch corresponding to the Lorentz-invariant structure 
"transverse" and the branch corresponding to the structure u>j° + k7 — "longitudi- 
nal" . These terms may be misleading in the fermion case because, in contrast to the 
plasmons with photon or gluon quantum numbers, these branches are not associated 
with transverse and longitudinal field polarizations. Hence the quotation marks. But 
better names were not invented, and using the words "transverse" and "longitudinal" 
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Figure 7: Plasmino spectrum in one loop. 

still makes a certain sense because the physical properties of "transverse" and "longi- 
tudinal" fermion branches are rather analogous to the physical properties of transverse 
and longitudinal gluon branches. 

The pattern of the quark spectrum is shown in Fig. 7. It is similar to the gluon 
spectrum with one important distinction — at k ~ 0, uj±(k) and uj\\(k) involve linear 
terms of opposite sign: 

k 

uJ q ± (k « ffT) = ^0 + - 

u «(k <^gT)=UJ Q -^ (5.9) 

where ujq is the plasmino frequency at zero momentum : 

a 2 T 2 

W ° 2 = 8 ° F (5 - 10) 
Thus, wjj(/c) first goes down and reaches minimum at some k*. The group velocity of 
the longitudinal plasmino at this point is zero. 

5.2 Landau Damping. 

The quoted one-loop results for the dispersive laws of transverse plasmons and plasmi- 
nos are gauge-invariant and stable with respect to higher-order corrections. The latter 
is not true, however, for longitudinal excitation branches 38, 1C]. We have seen that 



ujl loop (k) tends to the line uj = k exponentially fast at k 3> gT. This can be easily 
seen from the analysis of the dispersive equations for longitudinal branches which in 
the limit k 3> gT have the form 

w ,| + k ~ ^— In (5.11) 

" k (jj\\ — k 
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At k ^> gT, the solution exists when the logarithm is large and lo» — k is exponentially 



small. The logarithmic factor in Eq.( 5.11 ) comes from the angular integral 

d9 

uj — k cos 9 



(5.12) 



which diverges at u = k. This collinear divergence appears due to masslessness of 
quarks and gluons in the loop depicted in Fig.l. But quarks and gluons in QGP are 
not massless — their dispersive law acquires the gap ~ gT due to temperature effects. 
An accurate calculation requires substituting in the loops the dressed propagators. As 



a result, the logarithmic divergence in the integral (5.12) is cut off and the logarithmic 
factor in Eq. fl5.ll ) is modified: 



ln zrh \ ln (u - ky + opt* (5 - 13) 

Dressing of propagators amounts to going beyond one-loop approximation. Strictly 
speaking, to be self-consistent one should also take into account one-loop thermal 
corrections to the vertices (this procedure is known as resummation of hard thermal 
loops |39|]), but in this particular case these corrections do not play an important role. 
What is important is the cutoff of the logarithmic collinear singularity due to effective 
temperature-induced masses. 

Substituting ( |5.13| ) in ( |5.11 ), we see that the new dispersive equation does not at 



all have solutions with real oo for large enough k. This fact can be given a natural 
physical explanation. When k is small compared to gT, there is no logarithmic factor 
in the dispersive equation, the modification (|5.13 ) is irrelevant, and the dispersive law 



of longitudinal modes does not deviate from the one-loop result. Then the logarithm 
appears, the modification ( |5.13 ) starts playing a role, and, at some k** ~ gT, the 



longitudinal dispersive curve crosses the line uj = k. At this point the longitudinal 
polarization operator acquires the imaginary part due to Landau damping. 

In usual plasma, Landau damping is the process when propagating electromagnetic 
waves are "absorbed" by the electrons moving in plasma. In the language of quantum 
field theory, it is a 2 — > 1 process 

7* + e^e (5.14) 

In real time technique, that corresponds to a contribution to the imaginary part of 
the polarization operator so that both internal electron lines in the loop are placed 
on mass shell. At T = the standard Cutkovsky rules imply positive energies of all 
particles in the direct channel, and the imaginary part appears only due to the decay 
7* — ► e + + e~ . At T ^ 0, Cutkovsky rules are modified and both signs for energy are 
admissible. Physically, that corresponds to the presence of real particles in the heat 
bath so that the process ( |5.14 ) may go. 
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Figure 8: Plasmon spectrum. 

Also in QGP imaginary parts of polarization operators may acquire contributions 
due to Landau damping. The corresponding processes are 

9* +g^g, g* + Q^q, q* + g^q, q* + q^g,--- (5.15) 

where g*,q* are plasmon and plasmino collective excitations and q,g are the excita- 
tions with characteristic momenta of order of temperature (in this kinematic region, 
the dispersive laws are roughly the same as for tree quarks and gluons and the star 
superscript is redundant). 

The kinematic condition for the scattering processes (5.15) to go is that the fre- 



quency of collective excitations g* and q* would be less than their momentum. We 
have seen that the condition uj < k is realized indeed for longitudinal plasmon and 
plasmino excitations starting from some k** ~ gT. At k > k** , the Landau damping 
switches on and the dispersion law acquires an imaginary part. The imaginary part 
rapidly grows and very soon becomes of order of the real part. From there on it makes 
no sense to talk about propagating longitudinal modes anymore. The situation is the 
same as in usual plasma where longitudinal modes also become overdamped at large 
enough momenta and disappear from the physical spectrum H. 

The true pattern of plasmon collective modes is shown in Fig. 8. A similar picture 
holds for plasminos. 

5.3 Observability. 

One-loop dispersive curves are gauge- invariant. However, the question whether these 
curves are physically observable is, again, highly non-trivial. It is easy to measure 
explicitly photon dispersive curves in usual plasma — to this end, one should study 
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the propagating classical electromagnetic waves, measure the electric charge density 
(say, by laser beams) as a function of time and spatial coordinates and determine 
thereby the frequency and the wave vector of the wave. 

But there is no such thing in Nature as classical gluon field due to confinement, and 
no classical device which would measure the color charge density exists. Even more 
obviously, quark fields (which have Grassmannian nature) cannot be treated classically. 
Hence, one cannot really measure the energy and momentum of propagating colored 
waves in a direct physical experiment. 

What one can measure are correlators of colorless (in the first place , electromag- 
netic) currents. Modification of dispersive laws affects these correlators and that can 
be observed. However, a colorless current always couples to a pair of colored particles 
and, as a result, physical correlators involve some integrals of quark and gluon Green's 
functions which are related to quark and gluon dispersive characteristics only in an in- 
direct way. Also, thermal modification of vertices is as important here as modification 
of the Green's functions. 

However, there is one special point on the dispersive curves which can in principle 
be directly measured in experiment. This is the point on the longitudinal plasmino 
curve where its frequency acquires a minimal value and the group velocity turns to zero. 
Consider the problem of emission of relatively soft e + e~ or /U + /i~ pairs by QGP. In a 
"thermos bottle" experiment, one should make sure that the size of your thermos bottle 
is much less than the lepton mean free path. Otherwise, the leptons are thermalized 
and their spectrum is just Planckian. But in heavy ion collisions experiments, QGP 
is produced in small volume, the condition -L c har "C Lf™ e path is satisfied, and the 
spectrum of emitted leptons (and photons) can provide a non-trivial information on 
dynamic characteristics of QGP. 

The spectrum of soft dileptons was calculated in [ 40 1 . This is one of very few phys- 
ical problems we know of where the hard thermal loop resummation technique |39| 
should be used (and was used) at full length. The spectrum feels the effects due to 
quark and gluon interactions in the region ~ Pi+i- ~ gT — the spectrum at 

larger energies and momenta is the same as for the gas of free quarks. One particular 
source of soft dileptons is the process q\ — ► qT, + . The probability of this process 
has a "spike" for the momentum of q]_ and q*^ coinciding with the momentum k* on 
the longitudinal plasmino dispersive curve with zero group velocity. There are just 
many plasminos at the vicinity of this point and the phase space factor provides a 
singularity at Ei+i- = wj_(A;*) — W||(fc*) in the spectrum. Another spike comes from 
the process when a longitudinal plasmino with momentum k* annihilates with a longi- 
tudinal antiplasmino with the opposite momentum to produce a lepton pair with the 
energy 2u>||(/c*). 

Unfortunately, in the soft region, the main contribution in the spectrum is due to 
cuts. In other words, the most relevant elementary kinetic processes are not q* — > 
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q* + l + l~ or q* + q* — > Z + /~, but rather <?* + 5* — > g* + 5* + etc. The spikes 
actually have finite width due to collisional damping of collective excitations (the issue 
to be discussed in the next section), and one can hope to see only a tiny resonance on 
a huge background. Still, such a resonance in the spectrum is an observable effect. 

6 Damping Mayhem and Transport Paradise. 

6.1 Direct decay. 

In the previous section, we discussed the Landau damping contribution to the imag- 
inary parts of polarization operators and, correspondingly, to imaginary parts of dis- 
persion laws. It comes from the kinematic region u < k and is physically related to 
absorbtion of ingoing excitations by thermal quanta like in ( |5.15| ). However, we did 
not say a word about the contribution of direct decay processes g* — > q + q etc. in the 
timelike kinematic region u > k. 

That was with a good reason. On the one-loop level, the contribution of decay 
processes in imaginary parts is nonzero and is of order ~ g 2 T. Unfortunately, it 
depends on the gauge and, in some gauges, has even the wrong sign corresponding 
not to damping of excitations but to instabilities pdfl . The point is that such one- 
loop calculation is unstable with respect to higher-order corrections. It is very clear 
physically — quarks and gluons in QGP cannot be treated as massless but acquire 
dynamical masses due to thermal effects. And the decay of a plasmon or plasmino into 
two other collective excitations is not kinematically allowed. The only exception is the 
process g* — > q* + q* which in principle may go if 

N f >9c F -c v = 6 (6.1) 

But there are at most three light flavors in real QGP and decay processes can be safely 
forgotten. 

6.2 Collisional damping. 

Still, damping is there even in the timelike region u > k due to collisions g* + q* — > 
9* + Q* > Q* + Q* ~> 9* + 9* e tc This is also the main source of damping of transverse 
electromagnetic waves in usual plasma [|j . A rough estimate for collisional damping in 
QGP can be done very simply. 

The meaning of damping is the inverse lifetime of excitations. We have 

C ~ {rufe)- 1 ~ na tot (6.2) 
where n ~ T 3 is the density of the medium and a tot for excitations which carry (color) 
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Figure 9: 



charge has a Coulomb form 

tot 4 f dpi g 2 , v 

- -9]^- T 2 (6-3) 

We took into account the fact that the power infrared divergence for the integral of 
Coulomb cross section is effectively cut off at p± ~ gT due to Debye screening 0. As 
a result, we obtain the estimate 

(Q>9 „ g 2 T (6.4) 

Note that this value for the damping is unusually large. It is much larger than, 
say, the damping of photons in ultrarelativistic e + e~ - plasma. The latter can also be 
estimated from the formula ( |6.2| ), but a tot is now not the Coulomb, but the Compton 
cross section. The integral has now the form § eT dp\/p\ ~ ln(l/e) and the estimate is 

C 7 ~ e 4 Tln(l/e) < e 2 T (6.5) 

The question arises whether the new anomalously large scale ~ g 2 T has a physical 
relevance 0. We will return to discussion of this point a bit later. 

An accurate calculation of the damping of fast moving (k 3> gT) quark and gluon 
excitations in QGP has been done in [M, |ll[ (see also p2[). Consider the graph in 
Fig. 9 for the quark polarization operator where the lines with blobs stand for quark 
and gluon propagator dressed by thermal loops. The imaginary part of the whole 
loop in Fig. 9 depends on the imaginary parts of internal propagators. The imaginary 

n and due to Landau damping effects — see more detailed discussion below. 

12 Do not confuse this scale with the magnetic scale which is also of order g 2 T. The former is related to 
kinetic properties of the system while the latter refers exclusively to static phenomena. 

13 It can be shown jicj] that when calculating the leading contribution in £ in the kinematic region k ^> gT 
which is under discussion now, vertex corrections can be disregarded. 
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Figure 10: A graph contributing to Im H(k). Crossed lines are put on mass shell. One of 
the quarks in the internal loop has a negative energy. 

part of the gluon propagator due to Landau damping turns out to be of paramount 
importance. Physically, this contribution just corresponds to the scattering processes 
g* + q* — > g* + q* and g* + g* — > g* + g* as can be easily inferred if spelling out the 
exact gluon propagator as in Fig. 10 (there is also a similar graph with internal gluon 
loop). 

Dressing the quark propagator (indicated by the blob in Fig. 10) is important. If 
not taking it into account, the imaginary part of the propagator is 5 - function and 
the result for Im E(fc) has the form 

Im £(fc) = -^-ajp^cpTj [ dp 5(p - p cos 9) [ dcos9 [ ^ (6.6) 

OZ7T J J -l J „4 i JL,,,4 Eo 

The integral has a power infrared behavior at p > gT, but this divergence is cut off 
due to Landau damping effects [the second term in the denominator in the RHS of 
Eq.(^T(])]. Still, the integral in (|6.6| ) diverges logarithmically at p -C gT and the result 
for the damping inferred from Eq,( |6.6| ) is infinite. 

The crucial observation is that the dressed quark propagator does not have singu- 
larities on the real po axis. A self-consistent account of the collisional damping for the 
quark Green's function moves its singularities in the complex plane. As a result, 5 - 
function in the integrand is replaced by a smooth distribution with the width of order 
C ~ g 2 T. This smoothing cuts off the logarithmic singularity in ( |6.6[ ) at p ~ g 2 T. The 
other source for the cutoff could be provided by magnetic screening effects , but the 
latter is an essentially non-abelian phenomenon whereas the cutoff due to smearing out 
the 5 - function is a universal effect which occurs also in an abelian theory. 

To find the dispersion law, one should add Eq.( |6.6D to the one-loop result for 
and solve the dispersive equation. The solution is complex ui pole (k) = lo 1 loop (k) — i((k) 
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Figure 11: Polarization operator of soft plasmino. 



and the final result for 3> gT) is very simple. 

C q = a s c F T ln(CVfl) 
(° = a s c v Tln(Cy g ) 



(6.7) 



Only the coefficient of logarithm can be calculated. The constants C q , C 9 under the 
logarithm cannot determined. Actually, we will see shortly that these constants are 
gauge-dependent and cannot be defined in a reasonable way. 

Damping of excitations in another kinematic region k = ("standing" plasmons 



and plasminos) was studied in [39|. Here the vertex corrections are as important as 
corrections to the propagators, however an accurate analysis of |^]] shows that it suffices 
to take into account only one (hard thermal) loop corrections both in polarization 
operators and vertices. Consider for definiteness the damping of standing plasminos. 
The graphs contributing to the soft plasmino polarization operator are shown in Fig.ll. 
Using Keldysh technique one can derive in the limit of soft external momentum 



T> R (k) = -2ig 



d 4 qT 



It 

2 ^ 



2vr 4 q 

r R (k,q + k;q)G R (q + k)T R (q + k,k- 



(k;q)+ 



(6.8) 



14 A refined analysis which takes into account the modification of Im S(fc) for complex ui when one starts to 
move from the real w axis towards a singularity and which is beyond the scope of this lecture shows that the 
dispersive equation has actually no solutions and the singularity is no longer a pole, but a branching point 
p3 , 12, 44j. But this branching point is located at the same distance from the real axis as the would-be pole 
and brings about the same damping behavior of the gluon retarded Green's function G R (t) ~ exp{— £i} at 
large real times. There is a recent claim jl5) that in abelian theory G R (ui) does not involve singularities at 
all at finite distance from the real axis. This would result in a non-exponential decay of G R (t) at large real 
time — G R (t) ~ exp{— aTt ln(eTi)} (cf. Eq.(4.3). We do not think, however, that it is correct — such a 
form of G R (t) does not conform with a smooth behavior of G R (u>) on the real uj axis where our calculation 
of S(fc) is well under control and the cutoff due to a finite fermion width should be taken into account. 
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(G R and D R are the retarded Green's functions and T R are the retarded vertices). 
Subsituting here transverse and longitudinal parts of the gluon Green's function and 
one-loop vertices Y^ v and T M , calculating the integral (which in this case can be done 
only numerically), and solving the dispersive equation, one arrives at the result 

C«(fc = 0) = 1.43c F ^ (6.9) 
However, the gluon Green's function involves also a gauge-dependent part 

D R( - a Hq) = (a- 1) ^ (6 10) 

where a is a gauge parameter and an infinitesimal ie is introduced to provide for the 
right analytical properties. At first sight, this gauge-dependent piece should not affect 
the position of the pole. Really, one can use the Ward identities (which hold also at 
finite temperature for retarded propagators and vertices 46, |TI[) to derive 



S fi(tt) (fc) = -ig 2 Tc F (a - l)[G R (k)]~ 



Al 



2ir*q \l(qo+ie) 2 -q 2 } [(q - ief - q' 

G R (k + q)- G R (k)] [G R (k)]^6.ll) 



We see the presence of factors [G^(A;)] _1 both on the right and on the left. G R (k) is 
singular at the pole and [G R (k)]~ 1 is zero. One might infer from this that a gauge- 
dependent contribution to S^(pole) and hence to the corresponding solution of the 
dispersive equation determining the pole position is also zero. 

However, this is wrong [^7], 12]. The point is that the integral in ( |6.11 ) involves a 



severe power infrared divergence and is infinite at the pole. We have a thereby Ox co 
uncertainty. This uncertainty can be resolved by choosing k not exactly at the pole but 
slightly off mass shell. Then ^(k)]' 1 is not exactly zero and also the divergence 
in the integral is cut off by an off-mass-shellness. When the distance from the mass 
shell is small, the final result for Y, R ( a \k — ► pole) does not depent on this distance and 
is just finite. In the soft momenta region 



^) = - i(a ~^f TCF -r° 



This brings about a gauge-dependent part ~ g 2 T in the damping of soft plasminos. A 
similar analysis with the same conslusion can be carried out for plasmons. 

The same gauge-dependence shows up in the damping of energetic plasmons and 
plasminos, but in the latter case, this gauge-dependence is parametrically overwhelmed 
by the leading gauge-independent contribution (|6.7p involving the factor ln(l/<7). 



15 To determine the exponential asymptotics of G R (t) at large real times, we have to stay on the real co 



axis i.e. off mass shell [n2| 
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6.3 Observability. 

The observed gauge-dependence of the damping obviously indicates that it is not a 
physical quantity. This is definitely true at least for soft plasmons and plasminos 
where the gauge-dependent part and the gauge-independent part fl6.9| ) are of the same 
order ~ g 2 T. 

Indeed, it is not possible to contemplate a physical experiment where this quantity 
could be measured. That should be confronted with the case of abelian plasma where 
damping of electromagnetic waves is a perfectly physical quantity and can be directly 
observed by measuring the attenuation of the amplitude of a classical wave with time. 
But as we already noted, no classical gluon or quark waves exist. This observation 
refers also for damping of electron and positron collective excitations in the ultrarela- 
tivistic abelian plasma. It also has an anomalously large scale ~ e 2 T (with the extra 
logarithmic factor ln(l/e) for k 3> eT) and it also cannot be directly measured. 

One could try to observe the effects due to damping in gauge-invariant quantites 
like the polarization operator of electromagnetic currents. An accurate analysis which 
goes beyond the conventional hard thermal loop resummation technique and effectively 
resums a set of ladder graphs shows, however, that a self-consistent account of the 
corrections due to damping in the quark Green's functions and in the vertices results 
in the exact cancellation of the anomalously large scale ~ g 2 T in the final answer |pj| 1. 
(for a similar analysis with a similar conclusion in scalar QED sec [48]). 

However, there is a physical problem where the scale ~ g 2 T can in principle show 
up. This is the already discussed problem of lepton pair production in QGP. We 
have seen that the spectrum of leptons pairs invoves spikes associated with a special 
point with zero group velocity on the longitudinal plasmino dispersive curve. Going 
beyond the hard thermal loop aproximation and taking into account the effects due to 
collisional damping in the Green's functions and in the vertices would bring about a 
finite width for these spikes of order g 2 T and there is a principle possibility to measure 
this width. This problem has not been studied, and it is not clear by now whether the 
width of the spike can by calculated analytically and whether one can single out this 
spike out of the background. 

What one can say quite definitely is that this width crucially depends on modi- 
fication of vertices due to collisional effects and has nothing to do with the (gauge- 



dependent) position of the pole (or whatever the real singularity is |43|, 12, |44J]) of the 
quark and gluon Green's function. 

Thus, we are convinced that the latter is not a physical quantity probably even 
for energetic plasmons in spite of the fact that the leading contribution Q6.7|) is gauge- 



independent there. We just do not know how on Earth this quantity could be measured. 
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6.4 Transport Phenomena. 

There is a lot of kinetic phenomena in QGP which are physical and measurable. Indeed, 
nothing in principle prevents measuring the electric resistance of a vessel with QGP 
or studying the flow of QGP through narrow tubes. They do not depend, however, on 
the anomalous damping scale g 2 T, but rather on a much smaller scale 

(Tchar)- 1 ~5 4 Tln(l/ 5 ) (6.13) 

This scale already appeared in ( |6.5| ) determining the damping of electromagnetic waves 
in e + e~ - plasma. And it is also the scale which determines the mentioned physical 
effects of viscousity and electric conductivity, and many others — heat conductivity, 
energy losses of a heavy particle moving through plasma, etc. 

The appearance of the scale g A T\n{l/g) has a clear physical origin. All the men- 
tioned effects are inherently related to the rate of relaxation of the system to thermal 
equilibrium. The latter can be estimated as 

{r rel )- 1 ~ na trans (6.14) 



It looks the same as the estimate for lifetime ( |6.2| ) but with an essential difference — in 
contrast to ( |6.2[) , the estimate ( |6.14j ) involves the transport rather than the total cross 
section. The transport cross-section is defined as 



a trans = J ^(1- C OS0) (6.15) 

where 9 is the scattering angle. The factor (1 — cos 9) takes care of the fact that small- 
angle scattering though contributes to the total cross section, does not essentially affect 
the distribution functions n g (p), n q (p) and is not effective in relaxation processes. For 
the Coulomb scattering in ultrarelativistic plasma, the transport cross section is 

* trans ~g 4 [ *jLi~£ln± (6.16) 
J(gT) 2 p\ T z T z g 



Multiplying it by n ~ T 3 and substituting it in fl6.140 , the estimate ( |6.13 ) is reproduced 



Viscousity and all other similar quantities can be calculated analytically in the 
leading order (probably, magnetic infrared divergences prevent an analytic evaluation 
of these quantities in next orders in g, but this question is not yet well studied). It 
is interesting that Feynman diagram technique proves to be technically unconvenient 
here , and the good old Boltzmann kinetic equation is the tool people usually use (see 
e.g. @). 

Let us make two illustrative estimates which make clear how the relaxation scale 
( 6.13 ) depending on the transport cross-section ( |6.16|) arises. 



First, let us estimate the electric conductivity of QGP (it is the quite conventional 
conductivity, not the "color conductivity" which is sometimes discussed in the litera- 
ture, depends on the anomalous damping scale g 2 T, and is not a physically observable 
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quantity — we do not have batteries with color charge at our disposal). Suppose at 
t = the system was at thermal equilibrium so that the quark distribution functions 
are no(p,x, t = 0) = njr(|p|). When we switch on the electric field, the distribution 
function starts to evolve according to the kinetic equation 

<9n(p,x,t) <9n(p,x,i) dn(p,x,t) . , 

at + v 5x = eE dp + ' • • (fU7) 

Dots in RHS of Eq.( 6,17| ) stand for the collision term which becomes relevant at t ~ 
x ~ T re i ~ Tf ree p at h ~ [g 4 T ln(l / g)]^ 1 ■ Thus, the electric field brings about distortions 
of the distribution function which grow up to the characteristic value 

v eEvT 

dn ~ eE— noT free path ~ 41 ri , s 

T g*ln(l/g) 

At this point, collisional effects stop the growth (a particle drifting in external electric 
field collides with a particle in the medium, forgets what happened before, and starts 
drifting anew). The density of electric current in the medium is 

j = e/vi„d^~E iI g M (6.18) 

The coefficient between j and E gives the conductivity. 

Let us estimate now the energy losses of a heavy energetic quark in QGP. Of course, 
free quarks do not exist, but a physical experimental setup would be sending into the 
bottle with QGP a heavy meson Qq with open beauty or top. In QGP, the meson 
dissociates, and a naked heavy quark propagates losing its energy due to interaction 
with the medium. It goes out then on the other side of the bottle dressed again with 
light quarks, but not necessarily in the same way as before. When Mq 3> Kqcd-, this 
dressing does not essentially affect its energy. A heavy particle containing Q can be 
detected and its energy can be measured. 

Suppose a heavy quark is ultrarelativistic, but its energy is not high enough for the 
Cerenkov radiation processes to be important. Then the energy would be lost mainly 
due to individual incoherent scatterings. The mean energy loss in each scattering 
is AE ~ T (T — is a characteristic energy of the particles in heat bath on which 
our heavy quark scatters). The mean time interval between scatterings is Tf ree path ~ 
[g A T\n(l/g)]- 1 . We obtain 

-^~^^ = Ca 2 s THn(l/g) (6.19) 

UiX Tfree path 

This estimate turns out to be correct up to the argument of the logarithm which in 



reality is energy-dependent 5Cj. The numerical coefficient C = Attcf /3 was determined 
by Bjorken |5l[ 
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